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Abstract 

In this paper we analyze the spinning motion of the hovering mag- 
netic top. We have observed that its motion looks different from that 
of a classical top. A classical top rotates about its own axis which pre- 
cesses around a vertical fixed external axis. The hovering magnetic 
top, on the other hand, has its axis slightly tilted and moves rigidly 
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as a whole about the vertical axis. We call this motion synchronous, 
because in a stroboscopic experiment we see that a point at the rim 
of the top moves synchronously with the top axis. 

We show that the synchronous motion may be attributed to a 
small deviation of the magnetic moment from the symmetry axis of 
the top. We show that as a consequence, the minimum angular ve- 
locity required for stability is given by y/ A/jlHIx / Jf for I 3 > Ii and 
by yj JbHj (I 3 — I\) for I 3 < I x . Here, I 3 and I x are the principal and 
secondary moments of inertia, [i is the magnetic moment, and H is the 
magnetic field. For comparison, the minimum angular for a classical 
top is given by y/&nHIi/I% both for I 3 < I x and for J 3 > I x . 

We also give experimental results that were taken with a top whose 
moment of inertia I x can be changed. These results show very good 
agreement with our calculations. 

1 Introduction. 

1.1 The hovering magnetic top. 

The hovering magnetic top is an ingenious device that hovers in mid-air while 
spinning. It is marketed as a kit in the U.S.A. and Europe under the trade 
name Levitron™ [fll |2| and in Japan under the trade name U-CAS||. The 
whole kit consists of three main parts: A magnetized top which weighs about 
18gr, a thin (lifting) plastic plate and a magnetized square base plate (base). 
To operate the top one should set it spinning on the plastic plate that covers 
the base. The plastic plate is then raised slowly with the top until a point is 
reached in which the top leaves the plate and spins in mid-air above the base 
for about 2min. The hovering height of the top is approximately 3.0 cm above 
the surface of the base whose dimensions are about 10cmxl0cmx2cm. The 
kit comes with extra brass and plastic fine tuning weights as the apparatus 
is very sensitive to the weight of the top. It also comes with two wedges to 
balance the base horizontally. 

The hovering magnetic top serves as an excellent model to study many 
branches of physics. These include magnetism, electricity, hydrodynamics, 
mechanics of rigid body, wave mechanics, dynamical stability of non linear 
as well as linear systems, life time of quasi-stationary states and solid state 
magnetism to name just a few. For example, in the field of neutral parti- 
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cle trapping the hovering magnetic top is a vivid macroscopic device which 
demonstrates how it is possible to trap a neutral particle with spin in a static 
magnetic field. In fact, one of the recent examples of such traps was devised 
and applied to trap Rb 87 atoms and to form Bose-Einstein condensation 0, [J. 
The trap used in this experiment is based on a rotating magnetic field and 
was given the name TOP (Time-averaged Orbiting Potential) trap. Though 
the TOP trap consists of a time dependent field it nevertheless operates on 
roughly the same principles as the trap of the hovering magnetic top. In 
addition to understanding magnetic traps there are many other issues that 
can be addressed by studying the hovering magnetic top. Questions such as: 
How high can it hover?, How long does it hover?, What are the tolerances 
on the mass of the top and the tilt angle of the base? How the trap works 
in the quantum regime? and other questions have been studied recently 
[0, [7], ^], [| [TTJ, |TT[ whereas other questions such as how does friction affect 
the stability of the top are currently under study. 



1.2 The synchronous motion. 

The physical principles underlying the dynamical stability of the hovering 
magnetic top rely on the so-called 'adiabatic approximation' and have re- 
cently been discussed in several papers [|T2| , |I3| , |15|]. In these articles the 
top was modeled as an axially symmetric shape with its magnetization taken 
as a dipole pointing along the symmetry axis and situated at the center of 
mass. Observing the top as it hovers forced us however, to augment this pic- 
ture by assuming that the dipole is canted by a small angle A with respect to 
the symmetry axis of the top. We were motivated by the fact that the motion 
of the magnetic top as it hovers looked different from that of a classical top 
(but also by recalling that it would be forced by manufacturing tolerances). 
As a classical top slows down it starts to deviate from the vertical while it 
precesses with a frequency which is not synchronous to the frequency of the 
spin of the top around its axis. The axis of the hovering magnetic top is 
also canted to the vertical but the motion is synchronous. There is only one 
frequency involved. This can even be seen with the naked eye. The axis of 
the magnetic top rotates around the vertical but synchronously with the spin 
of the top around itself. An analogy, though not literal, is helpful here: The 
motion of the hovering magnetic top is reminiscent of the motion of the moon 
around the earth but that of a classical top is more like the rotation of the 
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earth around the sun. The astute observer watching the hovering magnetic 
top will note that the canting increases as the top slows down, even with 
the naked eye. It was more conspicuous to us when we observed it with a 
stroboscope preferably strobed at twice or three times the synchronous fre- 
quency. Using the latter frequency enabled also to measure the phase of the 
tilt. It is this way that we found this phase is correlated to the phase of A- 
the canting of the magnetic moment with respect to the axis of the top. We 
reach the same conclusion by observing in slow motion videos of hovering 
tops in which we changed A artificially. As will be shown in the following 
the canting A explains all the above observations. 

The present paper is in fact a study of the effect of this small canting 
A on the hovering of the magnetic top. In addition to explaining the above 
qualitative features we also derive quantitative results. In particular we show 
that the minimum angular velocity required for stability 0J m i n is given by 



UJr. 
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Here, ^3 and I\ are the principal and secondary moment of inertia, \i is the 
magnetic moment of the top and H is the magnetic field. This result is to 
be contrasted with the minimum speed of a classical top [[IB] which is given 



by ^JAnHIi/Il for both J 3 < l x and h> h- 

On the experimental side, we have measured the minimum speed of a 
hovering top whose moment of inertia I\ can be changed so as to cover values 
both below and above I 3 . These results show excellent agreement with our 
calculations. 



1.3 The structure of this paper. 

The structure of this paper is as follows: In Sec.(H) we carry out a detailed 
analysis of the synchronous motion. We first describe precisely what is meant 
by the synchronous motion and define our notations in Sec. (|2.1|) . In Sec.( |2.2| ) 
we study the conditions that must be met in order for the system to reach 
equilibrium. Next, in Sec.( |2.3| ) we carry out a dynamical stability analysis 
to find under what conditions the equilibrium solution found previously is 
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indeed stable. Finally, in Sec.( [2.4| ) we combine all the conditions that were 
found both from equilibrium considerations and dynamical stability consid- 
erations, and arrive at the expressions for the minimum angular velocity. The 
description of our experiment and a comparison of the results with the the- 
oretical calculations from the previous section is given in Sec.(|3|). In Sec.(f|), 
we summarize our results and discuss possible extensions of the calculations 
presented in this paper and other related subjects. 



2 Analysis of the Synchronous motion. 
2.1 Description of the Problem. 

As the top hovers above the magnetized plate, the magnetic lift force ex- 
erted by the base is balanced by the top weight. Therefore, to simplify 
further calculations we assume that gravity is zero and the magnetic field 
is homogeneous. Further, we disregard the translational motion of the top 
and assume that it move around a fixed point. This point is taken as the 
common origin of the moving and fixed system of coordinates. With these 
simplifications the synchronous motion can be modeled as shown in Fig. ([!]). 

The figure describes a top that rotates rigidly around the vertical z axis 
with angular velocity uj. The principal axis of the top (axis ii) makes an 
angle 8q with the vertical axis. The whole system is in a uniform magnetic 
field H, pointing downward. The top possesses a magnetic moment dipole, 
/i, that makes an angle A with its principal axis. 



2.2 Equilibrium considerations. 



We first determine the angle between the top axis h and the vertical axis 
z from torque equilibrium. The angular momentum of a symmetric top has 
the formjlTj 

t - dn /1X 
L = I 3 u 3 n + I 1 nx— (!) 
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where us is the component of the vector angular velocity cD in the direction 
of the n axis. Since h is fixed in the body, 



dh A . . 

- = cxn, (2) 

and therefore 

L = (J 3 - fx) ^ 3 n + JiW. (3) 

For the synchronous motion, the gyroscopic torque in the co-moving frame 
must compensate the magnetic torque, hence 

uj x L = n x H. (4) 

Thus, /x has to be coplanar with z and n, and one obtain, with = u cos #o, 
the equilibrium condition 

(J 3 - Ji)u; 2 sin 6^0 cos 6^0 = fiH sm(6 + A). (5) 

We assume that the angle A is small. Then, during its motion the top 
axis n stays close to the z axis, and we may use small angle approximation 
in Eq. fl5|) and express 8q in terms of u> as follows: 

- - " ffA (6) 



(/„ - h)u> 2 - liH 

In order to understand what Eq. (Bp means we must resort to a more realistic 
description of the hovering top. If we take into account air friction and the 
inhomogeneity of the magnetic field, the following picture emerges: As the 
top hovers, it experience air friction which causes uj to decrease with time. 
Thus, according to Eq.@ 6 increases provided that J 3 > I\ (a disk-like top). 
This, on the other hand, decreases the repulsive force on the top which is 
caused by the field gradient. Thus, the top rebalances itself at a lower height. 
This process continues until the field gradient reaches its maximum value, 
that occurs at a height equal roughly to half the size of the base. At this 
point the magnetic field can no longer support the top and the latter falls 
down vertically. In our experiments we could clearly see that 6$ increases 
gradually by the naked eye. We conclude that when 9 is 'large enough' 
the top will fall down. We denote by 6q the maximum allowed angle. The 
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angular velocity which correspond to this angle is, according to Eq.(|6]), given 
by 




To estimate 9q we recall that the tolernace on the weight is equal to the 
toleance allowed on the lift. It follows that[|TU 



dm (6*Y 



m 2 

Since the tolerance on the weight is about 2% (see Ref.|lO|) it follows that 
9q ~ O.lrad. Although we did not measure A directly, our experiments, 
including those in which we have artificially changed A indicate that A/#q <C 
1. It is therefore justified to neglect this term in Eq.(|7]) and we arrive to an 
approximate expression for the minimum speed when J 3 > namely: 



For a rod-like top (J3 < I\) the synchronous motion is also possible, but 
now A must be more negative than — 8q. This is because in this case the 
direction of the centrifugal torque is reversed, and for equilibrium to exist the 
magnetic torque must also be reversed. The condition that A < — 9 gives 
uj > which signifies that as far as the equilibrium condition is concerned 
there is no minimum speed for this case. Following similar arguments as 
before it can be easily deduced that in this case 9q decreases with time and 
the top slowly gains height as it slows down due to air friction. This is 
opposite of what we have found for a disk-like top. 

Note also that the rod-like top exhibits diamagneic-like behavior. Namely, 
the differential susceptibility is in this case negative as can be seen from 
Eq.©. 

Summarizing our results so far we conclude that as far as the equilibrium 
condition is concerned the minimum allowed angular velocity depends on the 
eccentricity of the top and is given by 



stat _ J J for a disk-like top (7 3 > h) ^ 



for a rod-like top (J 3 < 7i), 
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and is valid both for clockwise and counter-clockwise directions of the spin. 
Note that Eq. (§) is not the all answer because the dynamical stability should 
also be considered. This is done in the following section. 

2.3 Dynamic stability considerations. 

2.3.1 Overview 

Up to this point we have studied conditions that were derived directly from 
the equilibrium condition. In this section we study the conditions that must 
be met in order to sustain the synchronous motion. We start by reformulating 
the equations of motion in terms of the Lagrangian formalism and Euler's 
angles, and then add to the synchronous motion a small perturbation. We 
then arrive at a set of equations for the perturbational part and find the 
conditions that are required for the perturbation to be oscillatory which is 
the condition that the unperturbed motion is stable. 

2.3.2 The Lagrangian and the equations of motion. 

We denote the coordinate system fixed in space by (x, y, z) and the other 
that is fixed in the body by (1,2, n), as depicted in Fig.(0). The magnetic 
moment /x makes an angle A with the n axis and lies in the 2-n plane. Thus 

= ft (n cos A — 2 sin A^ . 

The magnetic field is uniform and points along the —z direction, i.e. 

H = —Hz = —H (1 sin ifj sin 9 + 2 cos ip sin 9 + n cos 9^j , 

where 0, 9 and ip are the first, second and third Euler's angles, defined in 
Fig.(0). Therefore, the magnetic energy for this configuration has the form 

E mag = — /x ■ H = fiH (cos 9 cos A — sin A cos ifj sin 9) . 

The kinetic energy, expressed using Euler's angle, is 

E km = \h (V> + cos 9) 2 + (9 2 + 2 sin 2 O) . 
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The Lagrangian L is given by the difference between the kinetic energy Eki n 
and the magnetic energy E mag , 

L = h 3 (V> + 4> cos 9) 2 + [O 2 + <j> 2 sin 2 6) (g) 

—fj,H (cos cos A — sin A cos ip sin 0) . 

The dynamics of motion for this system is then governed by the Lagrange 
equations of motion 

d fdL\ _ dL _ 
di Vd0 j ~ ~d0 ~ 

d_ fdL\ _ dL_ _ 
~dt \di) ~ !ty ~ 
With the Lagrangian given in Eq.(|9]) one obtains 

1\0 + (I 3 — Ii)(p cos + I^i/j (J) srn 9 = fiH (sin cos A + sin A cos ip cos 0) 
l^cf)0 sin — 7 3 yip + (j) cos 0^j = fiH sin A sin ip sin 
(la cos 2 + 1 1 sin 2 0) (j) + I 3 ip cos = L z = Const. 

(11) 

2.3.3 The Stationary Solution. 

A stationary solution to the set Eqs. (|TTD is found by setting 

0(t) = O = Const. (12) 

ip(t) = ipo = Const. 

Substituting these into the third equation of the set Eqs.([TT|) and integrating 
over time yields 

^{t) =tu = Const. (13) 
Using the second equation we find that 

V>o = 0, (14) 
while the remaining (first) equation gives 

(I 3 - iiV 2 sin 6» cos O = fiHsm(0 o + A). (15) 
Note that Eq.(pT5|) is identical to Eq.(|5|), which was derived in Sec.(E72|). 
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2.3.4 Perturbing the stationary solution. 

We now perturb the stationary solution by adding small variations, 

4>(t) = 4> s (t) + 5<p(t) = u + 5<p(t) 

ip(t) = ^ s (t) + <fy(t) = + <fy(t) (16) 
9(t) = 3 (t) + 60(t) = 9 + S9(t) 

Substituting this these into Eqs.fllTD and keeping only first-order terms yields 

hse + (h - h)uj5<j)sm(2e Q ) 

+(J 3 - h)u 2 59 cos(2# ) + hSipu sin 6 = fiH59 cos(A + 9 ) 
I 3 (5tp + 5(f) cos 9 ) - I 3 59u sin 9 + fiH sin Adtp sin 9 = (17) 
— 2u>(I 3 — h)59 sin 9q cos 9q + I 3 8ip cos 9q 
+50 (/ 3 cos 2 # + h sin 2 O ) = 

Using Eq.(|Uj) we express A in terms of 9 and to simplify further calcu- 
lations we now restrict ourselves to the case where 9q is a small angle and 
work to second-order in 9q. Under the small-angle approximation Eqs.(|17|) 
becomes 

hS9 + 2(/ 3 - h)u6<p9 + (I 3 - h)u 2 69 (1 - 29 2 ) 
+I 3 5ipuj9o = /^H59 (1 - 0.5{3 2 9%) 

i 3 5i> + i 3 S4> (i - ^ 2 ) - hseue 

+(oo 2 (I 3 -I 1 )- f iH)9 2 6^ = 
-2u9 Q {I 3 -I x )59 + I 3 8$(l-\9t) 
+I 3 64>(l-9 2 ) + I 1 64>9 2 = 0. 

where (3 = (J 3 - h) uj 2 /fiH - 2. 

Since Eqs.([l8|) are linear in the perturbations and homogenous, the gen- 
eral solution of these equations is a linear combination of exponential func- 
tions. We therefore set 

89{t) = 59 e xt 
S(j)(t) = 5(p e xt 
5if){t) = 5ip e xt 

and substitute these into Eqs. (|l8|) . The result is the matrix equation for the 
determination of the eigenvalues A 



10 





' se - 




' " 


A 



















where A is the matrix given by: 



\ 2 h - (j,H{l - 0.5(3 2 9 2 ) 
+ (/ 3 - h)u\l - 29 2 ) 

-\I 3 w9 



2X(I 3 - h)u9 

A 2 / 3 (i - \el) 

< A 2 / 3 (l - £ 2 ) > 
V +A 2 /i^o 2 J 



\I 3 ^9 
\ 2 h 



+ [u 2 (J 3 - h) - fiH] 9 
o 2 ) 

(19) 



A 2 / 3 (l - ^ 2 ) 



For a nontrivial solution to exist, the determinant of the matrix must vanish. 
After a few algebraic steps we find 



det A = -\ 2 h9l (a\ A + b\ 2 + c) + O (^) = 0, 



where 



a = I 2 

b = -2h[nH-u 2 {h-h) 
c = [fiH-uj 2 (h-h) 12 



- oo 2 {h-2h){2h-h) 



Thus, we find a doubly degenerate zero eigenvalue corresponding to the cyclic 
coordinate 0, which enter the equations of motion only through its time 
derivative. The other eigenvalues are the roots of a quadratic equation in A 2 . 

For the perturbation to be bounded for all times, the non-vanishing eigen- 
values must be purely imaginary. Thus, the two A 2 roots must be both real 
and negative. The condition for these roots to be real is that b 2 — Aac > 0. 
This results in the inequality 



»H 1 fh 

uj 2 h 4 Ui 



(20) 
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which is the well-known stability condition of a classical top with one point 
fixed flT6|. For both roots to be negative, their product must be positive (i.e. 
c/a > 0) and their sum must be negative (—b/a < 0). Since both a and c are 
squares of real numbers, their ratio is positive. We thus find as the second 
condition 



yH ^ 1 // 3 \ 2 (I 



It may be easily verified that for a physical top, for which 73 < 2ii, the 
inequality Eq. (|2lD is satisfied whenever Eq.(pO|) holds. Therefore, the overall 
condition for the stability of the synchronous motion is given by Eq.(|20|) 
alone, namely 

dyn _ 



^ (22) 



both for I 3 > I\ and for J 3 < I\. 



2.4 Conclusions. 



Hitherto we have found separate conditions for the minimum angular speed, 
one pair given in Eq. @ arising from equilibrium considerations, the other one 
given in Eq.(^) originating from dynamical stability considerations. Since 
these conditions must be satisfied simultaneously, we now take the union 
of these conditions. Taking into account the fact that for a disk-like top 



fiH/(I 3 — Ii) is larger than ^jAfiHIi/ If, we find that the union of these 
conditions amounts to the following expressions for uj min : 



HH 



h 



T 2 



for a disk-like top (I3 > Ji) 
for a rod-like top (I3 < Ii). 



(23) 



According to this result, there should be an abrupt change in w m j n in the 
vicinity of a sphere-like top (I\ = 1%) even for a vanishingly small value of the 
angle A. This abrupt change is also found in our experiments, as is shown 
in the next section. 
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3 Experimental Results. 



To test our calculations we have built a magnetized top whose moment of 
inertia I\ can be changed so as to cover values both below and above 1%. 
We have measured the minimum angular speed of the top for several values 
of I1/I3. Fig. (|) shows (in bullets) the reciprocal of the measured angular 
speed normalized to the minimum angular speed of a pure disk top versus x = 
(I 3 — To compare these results to our calculation we rewrite Eq . fl23|) 

as an equation for the normalized reciprocal angular speed 0J min (1) /u) min (x) 
in terms of x. The result is 



UJ 7l 



/ 2x 
x + 1 
lx + 1 



< x < 1 
-1 < x < 0. 



Fig.(|3p shows uJrain (1) / 'u min (x) as a solid line. It is clear from the figure that 
the experimental results support our theory. 

The two main features of our theoretical results are: 1) the divergence of 
u) min for disk-like tops as they approach spherical shape and 2) the disconti- 
nuity in uJmin as one goes over to rod-like top. 1) Unfortunately, we could not 
go beyond u > 2.5u min (l) because the maximum spin frequency is limited 
by the dynamic stability of the top in the trap|14|. It is therefore impossible 



to observe experimentally the above divergence. 2) Though we have only two 
measured points at the rod-like part of the graph it is conspicuous that they 
verify the difference in the origin of the instability for disk- and rod-like tops 
that our theory shows. 



4 Discussion. 

In this paper we have analyzed one possible cause for the synchronous motion 
and attributed it to the small difference between the direction of the principal 
axis of inertia and the direction of the magnetic moment. We have shown that 
as a consequence, for a disk-like top the minimum angular velocity required 
for stability differs from that for a classical top (with one point fixed) and is 
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given by: 



Mr, 



for J 3 < I x 
for I 3 > h 



V 4 - A 

where J3 and I\ are the principal and secondary moment of inertia, \i is the 
magnetic moment of the top and H is the magnetic field. This result is found 
to be in a very good agreement with our experiments. 

We remark that in the calculations presented in this paper we have ne- 
glected some aspects which are relevant to the behavior of the top. For 
example, the dipole moment of the top carries with it a minute amount of 
intrinsic spin of the order of S ~ h\xj \xb where h is Planck's constant and 
Hb is the Bohr's magneton. The existence of the spin contributes an addi- 
tional term to the total angular momentum written in Eq.@. In this case 
the symmetry between clockwise (CW) and counter-clockwise (CCW) rota- 
tions breaks down, and one finds different values for the minimum angular 
speed for each direction. Due to the smallness of the spin with respect to the 
orbital angular momentum, the CW-CCW asymmetry becomes pronounced 
only when the size of the top is of the order of tenths of microns. Our cal- 
culations show that when the sense of rotation is roughly antiparallel to the 
direction of the spin the minimum angular speed increases, whereas when it 
is roughly parallel to the spin, the minimum angular speed decreases. This 
conclusion immediately raises the question whether it is possible to hover 
a top using its intrinsic spin alone without the need to spin it. To study 
this we assume the top to be endowed intrincily with spin proportional to 
its magnetic moment. We find that it may be possible to hover a top having 
only this spin with no additional angular moment um[|i~8| . 

Another interesting point is what happens if we relinquish the requirement 
that #0 is small. Note that the equilibrium equation Eq.([|) is that of the 
well-known Stoner and Wolfarth model (SW)fll9| [20fl , extensively used to 
described the hysteresis of single domain ferromagnetic particles. In that 
case at most only two of the equilibria are stable. We find, however, in our 
case the unstable SW becomes stable under some conditions, being stabilized 
by the dynamics. A word of warning is due here when applying these results 
to the hovering top. We have in this work considered only the rotational 
degrees of freedom, and what is stable under these is not necessarily stable 
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any more in an inhomogeneous field necessary to trap the top. We believe 
however that the small angle solution is still stable there under the same 
conditions that are given by us JT4J for A = 0. 
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Figure 1: The synchronous motion: The top revolves rigidly around the 
vertical z axis with angular velocity u. The principal axis of the top (axis 
n) makes an angle 6q with the vertical. The whole system is in a uniform 
magnetic field, H, pointing downward. The top possesses a magnetic moment 
dipole, fi, that makes an angle A with its principal axis. 
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Figure 2: Definition of Euler's angles. 
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Figure 3: Reciprocal of the normalized minimum angular velocity versus 
top shape. The bullets are the measured results and the solid line is the 
calculated results. 
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